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The tunneling spectrum of an electron and a hole in metal-superconductor-metal junctions is 
computed using the Blonder- Tinkham-Klapwijk theory. The incident and the outgoing currents 
finally balance each other by an interface charge inside the superconductor and metal junction. The 
present computation shows a more abundant structure compared to that on a metal-superconductor 
junction, such as the resonance at bias voltages above the energy gap of the superconductor. The 
density of the interface charge shows a quantum-like oscillation. 

PACS numbers: 



I. INTRODUCTION 



The Andreev reflection on metal-superconductor junc- 
tions has a long history and still attracts great interests of 
scientists recently Q~9< The tunneling spectrum of elec- 
trons on a metal-superconductor junction is sensitive to 
the energy gap of the superconductor, thus becoming an 
important technique to measure the gaps and gap proper- 
ties of superconductors. When an electron transmits into 
a superconductor to form a Cooper pair there appears a 
hole reflection, which is called the Andreev reflection. A 
commonly used description to the Andreev reflection is 
the Blonder- Tinkham-Klapwijk(BTK) theory 0, which 
take the interface in a metal -superconductor junction as 
a 6(x) potential barrier. Wei, Dong and Xing studied the 
tunneling spectrum in metal-superconductor-metal junc- 
tions (MSMJs) H using the BTK theory. They found 
unreasonable results that the incident and outgoing cur- 
rents in both sides do not balance each other. Thus they 
claimed that the BTK theory are not suitable for the 
Andreev reflection in MSMJs and treated them in the 
Landaucr-Buttiker formalism [5] . 



The authors of this paper think that since the BTK 
theory uses the standard quantum mechanics and did 
not do any unreasonable approximation it should be 
still working well in MSMJs. The unbalanced currents 
must cause a charge accumulation on the interface in the 
MSMJs. These interface charges change the potential on 
the metals and finally reaches a balance. This dynamic 
process can be simulated by solving the time-dependent 
Schrodinger equation. In the present paper a simpler 
model is set up to describe the final balanced state of 
the MSMJs. The differential conductance of MSMJs is 
computed in this model. 
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FIG. 1: A metal-superconductor-metal junction with an inci- 
dent electron and an incident hole from the left side and the 
right side respectively under a bias voltage. The metal leyer 
and the superconductor are labeled by M and S, respectively. 



II. FORMALISM 



The metal-superconductor-metal junctions are shown 
in Fig.l, where two thin insulating layers exist inside 
the two junctions becoming two potential barriers. An 
electron and a hole are incident into the junctions under 
a bias voltage from the left side and the right side, re- 
spectively. Due to the two barriers the electron and the 
hole are partially reflected and partially transmitted into 
the superconductor (SC). In the SC a pseudo-particle can 
only propagate above the energy gap A. When the bias 
voltage is smaller than A only Cooper pairs can propa- 
gate in the superconductor. In this case there will be a 
hole reflected on both sides, respectively, which are called 
Andreev reflection in literatures. 

A key idea of this paper is that the metal- 
superconductor interface at the right side will be charged 
due to the supercurrent in the SC. The tunneling process 
of an electron and a hole on the MSMJs is described by 
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FIG. 2: The energy dispersion on the MSMJs in the case 
of eV > A. The arrows label the propagating directions of 
quasi-electrons (red) and quasi-holes (blue). The charges on 
the interface are labeled by a column of which can also 
be positive. 



SC and propagate as a quasi-particle. It is then partially 
reflected by the following interface and partially tunnels 
into the following metal. So does an incident hole from 
the right side. The wave functions of the electron and 
hole are written as 
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the Bogoliubov-de Gennes (BdG) equation given by[ 
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where V(r) = eV, —eV in the left and right metals, re- 
spectively, and in the SC. On the interfaces V(r) = 
U5{x). eV = eV + eSV is the potential energy at 
the right side which is slightly changed by the interface 
charge. The continuity equation of the BdG equation is 
given by 



at 



V-(J Q (r)+J s (r)) = 



(3) 



where pq(y and 3q(t) = ^Im{u*\7u + v*Vv) are the 
charge density and the current density of quasi-particles, 
respectively. The supercurrent density is defined by 
V • Js(r) = -|/m(Au*(r)t)(r)). It can be proved 
that quasi-particle current is continuous along the whole 
MSMJs, but the supercurrent is obviously discontinuous 
since no supercurrent exists in a metal. This charges the 
interface like the charging process on a capacitor with a 
bias voltage. This can be seen by integrating the con- 
tinuity equation over a slap volume across the interface, 
resulting in 



da(t) 
dt 



= Js(x = L) 



(4) 



where a(t) is the surface charge density on the interface 
and Js{x = L) is the supercurrent at the interface be- 
tween the SC and metal II. The interface charge density 
reaches saturation when the current densities on both 
sides of the interface balance. The interface charge dis- 
tribution changes the potential energies of electrons and 
holes in metal II thus the outgoing current is changed 
and finally reaches balance with the incident current. 

As seen from Fig. 2 when the bias voltage satisfying 
eV > A an incident electron in metal I tunnels into the 



are the quasi-eletron and 



where ( U+ ) and ( U 

quasi-hole wave functions, respectively. They are given 
by 
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The wave vectors are determined by 
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In the case of eV < A, q± become complex numbers, 
so that the wave functions of electrons become damping 
traveling waves and Cooper pairs appear. 

The constants a, b, c, d, e, /, g, h in to ([JJ are ob- 
tained from the boundary conditions of wave functions 
on the interfaces, where the wave functions are continu- 
ous but the first derivatives of them have jumps due to 
the 5(x) barriers, i.e., 
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where Z 



2mU 



h2k and L is the thickness of the SC. Finally 
the current density, the super-current and the differential 
conductance are computed with these constants, which 
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FIG. 3: Charge current densities inside the SC without(left, 
a = 0.) and with(right,a = 0.24 x 10~ 4 e) interface charging. 
The parameters are set to fi = 0.5, A = 0.001/i, Z — 1, kpL = 
10000, eV = 1.8A. Circles label the incident and outgoing 
currents in the metals. 



are given by 
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FIG. 4: The incident and outgoing current densities versus 
different interface charge densities at eV — 1.8A. Other pa- 
rameters see Fig. 3. 
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FIG. 5: The normalized differential conductance without(left 
panel) and with (right panel) interface charges versus bias 
potential. Parameters see Fig.3. 



III. RESULTS OF COMPUTATION 



To view the effect of the interface charges we first com- 
pute the current densities in the MSMJs with and with- 
out interface charges. As shown in Fig.3(left), the charge 
currents is continuous but the incident and outgoing cur- 
rent do not balance with each other. It is natural to 
think that the interface at the right side will be charged 
until saturation so that both the currents balance. This 
is true. As seen in Fig.3(right) the two currents equal 
to each other after fully charged. It is interesting to 
note that the super-currents are always complementary 
to the charge currents in both cases so that the total 
currents are always constant in the SC. In addition, the 
quasi-particle current and the super-current with inter- 
face charges become symmetric about the center of the 
SC. This result convinces us that the interface charges 
have a true effect on the currents in the MSMJs. 

Further computation shows that the the incident 
and outgoing currents in fact oscillate on the interface 
charges. There exists a series of interface charge densi- 
ties which balance the incident and outgoing currents, as 
shown by the matching points in Fig. 4. The physics of 
this periodicity is still an open question. One may find a 
quantum behavior of the interface charge densities. 

The tunneling spectrum, i.e. the normalized differen- 
tial conductivities at different bias voltages without and 
with interface charges are shown in Fig. 5. Without in- 
terface charges the incident and outgoing currents match 
each other only under small bias voltages. After fully 
charged the two currents overlap in the whole bias region. 
Therefore, the conservation of currents in the MSMJs is 



realized by means of the interface charging. In addition, 
the probability flux \a\ 2 + \b\ 2 + \ f\ 2 + \g\ 2 = 2 is always 
conserved for different bias voltages. 

Finally the tunneling spectra for the differential elec- 
tric conductance (DEC) across different interface po- 
tential barriers are computed, as shown in Fig.6(lcft). 
Across a vanishing barrier the DEC approaches to a con- 
stant 4, which originates from the electron and hole in- 
jection. This is very different from that across a metal- 
superconductor junction, where the DEC drops gradually 
from a constant value to about half at the bias voltage of 
the energy gap(A). The main difference between these 
two cases is that there is a hole injection from the right 
side across the MSMJs. As the interface barrier increases 
the DEC shows a few sharper and sharper resonances. 
Their positions depends on the width of the supercon- 
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FIG. 6: The differential conductance (left) and the interface 
charge density (right) versus the bias voltage for different in- 
terface barriers. Other parameters see Fig.3. 
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ductor inside the MSMJs. Below the bias voltage eA the 
DEC decreases rapidly similar to the case across a metal- 
superconductor junction. This provides a measurement 
to the energy gap of the superconductor. The interface 
charge shows an interesting transition at the gap when 
the interface barrier reaches 4, see Fig.6(right). This may 
be taken as a signal that the resonance is fully set up. 

In summary, the tunneling spectrum of an electron and 
a hole on MSMJs is computed through the Bogoliubov- 
de Gennes (BdG) equation. It shows a more abundant 
structure compared to that on a metal-superconductor 
junction, such as the resonance above the gap voltage. 
An important effect occurring in the present device is 
that the incident current and the outgoing current bal- 



ance each other only when there is an interface charge in- 
side the superconductor and metal junction. The density 
of this interface charge shows a quantum-like oscillation, 
where a scries of densities balance the currents. 
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